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Summary 


A  set  of  finite  difference  calculations  were  performed  to  model  the  effects  of  a 
quarry  bench  on  the  excitation  of  seismic  signals  from  a  quarry.  This  work 
complement?  other  S-CUBED  research  on  the  nonlinear  behavior  of  the  quarry 
blast  source  (spall  from  the  quarry  face)  and  interference  effects  from  ripple-fired 
explosions.  We  find  that  the  excitation  of  seismic  phases  depend  strongly  on  the 
location  of  the  explosion  source  near  the  quarry  face.  In  particular,  we  find  that: 

1.  Seismic  signals  from  point  explosion  (dilatational)  sources  located  behind  the 
quarry  face  can  be  reduced  substantially  in  amplitude,  while  seismic  waves 
from  sources  on  the  quarry  floor  may  be  amplified.  The  source  amplitude 
variations  occur  for  sources  within  about  one  bench  height  of  the  quarry  face. 
The  numerical  simulations  presented  in  this  study  can  be  used  to  estimate 
and  correct  for  this  effect. 

2.  An  extended  source  in  the  standard  configuration  (charges  placed  a  distance 
of  1/2  the  bench  height  from  the  quarry  face  and  extending  to  the  depth  of 
the  quarry  floor)  averages  these  variations  leading  to  a  source  that  is  reduced 
slightly  in  amplitude  relative  to  a  source  in  a  half  space. 

3.  The  amplitude  variations  occur  because  a  dilatational  source  such  as  an 
explosion  is  strongly  affected  by  the  presence  of  a  free  surface  close  to  the 
source.  Seismic  signals  from  point  forces  are  affected  very  little  by  proximity 
to  the  quarry  face. 

4.  These  results  are  nearly  independent  of  the  takeoff  angle  of  the  seismic  signal. 

5.  Rayleigh  wave  amplitudes  are  reduced  by  up  to  40  percent. 


# 
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1.  Introduction 


In  our  recent  annual  report  (Barker,  et  al.  1993),  we  modeled  several 
important  aspects  of  the  quarry  blast  source  including  horizontal  and  vertical 
spall  from  the  quarry  blast  face  and  the  effect  of  superposition  of  ripple-fired 
explosions.  We  also  examined  the  effectiveness  of  a  large  quarry  blast  source  for 
hiding  a  clandestine  nuclear  test.  In  this  report,  we  examine  the  effect  of  the 
location  of  the  source  within  the  quarry.  Stevens,  et  al.  (1993),  found  that  surface 
waves  from  underground  nuclear  explosions  could  be  reduced  substantially  in 
amplitude  by  the  presence  of  a  vertical  boundary  close  to  the  explosion.  This 
brings  up  the  question  of  how  the  quarry  face  affects  the  seismic  signal  from  the 
quarry  blast,  whether  it  causes  any  of  the  seismic  signals  from  the  quarry  blast  to 
be  amplified  or  reduced  in  amplitude,  and  whether  it  leads  to  any  characteristics 
of  the  seismic  signals  that  could  be  used  to  identify  the  source  as  a  quarry  blast. 

To  answer  these  questions,  we  have  performed  a  series  of  finite  difference 
calculations  for  an  explosion  located  at  a  suite  of  locations  near  a  quarry  bench. 
We  also  perform  a  theoretical  analysis  and  find  an  analytic  solution  to  this 
problem  valid  for  wavelengths  much  larger  than  the  height  of  the  quarry  face, 
which  includes  nearly  all  seismic  signals  of  interest. 
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2.  Numerical  Simulations  of  Explosions  Near  a  Quarry  Bench 

2.1  Methodology 

An  explosive  charge  in  a  standard  configuration  for  a  quarry  blast  is  shown  in 
Figure  1.  The  seismic  signals  from  this  source  result  from  a  linear  component,  the 
signals  due  to  a  dilatational  source  in  this  structure;  and  a  nonlinear  component, 
the  nonlinear  motion  (spall  and  ballistic  motion)  of  the  material  broken  from  the 
quarry  face.  To  calculate  the  linear  component  to  the  quarry  blast  source,  we  use 
the  approach  described  by  McLaughlin  and  Jih  (1988),  using  the  elastic 
reciprocity  theorem  (e.g.  Aki  and  Richards,  1980)  to  find  the  far-field 
displacements  due  to  a  point  explosive  source  by  calculating  the  dilatation  at  the 
desired  source  point  caused  by  an  incident  plane  wave  (see  Figure  2). 

This  reciprocal  approach  has  several  advantages  over  the  conceptually  more 
direct  approach,  which  would  be  to  locate  explosive  and  force  sources  at  source 
locations  near  the  bench  and  compute  the  resulting  far-field  wavefields.  The  first 
advantage  is  that  the  reciprocal  approach  yields  the  far-field  wavefield.  In  order 
to  find  the  far-field  signals  from  the  direct  approach,  one  must  either  use  a 
numerical  grid  sufficiently  large  to  propagate  the  signals  several  wavelengths  to 
the  far  field  or  devise  a  continuation  scheme  similar  to  the  one  used  in  numerous 
S-CUBED  reports  to  find  the  far-field  waves  from  nonlinear  sources  (e.g.,  Barker 
and  Day,  1990).  For  5  Hz  signals  in  a  material  with  a  P  wave  velocity  of  5000 
m/sec,  the  wavelength  is  1000  meters,  so  the  grid  would  need  to  be  at  least 
several  km  in  each  direction  to  reach  the  far  field  and  several  more  km  to 
minimize  interference  from  boundary  reflections.  Adequate  description  of  the 
topographic  bench  requires  a  grid  spacing  less  than  5  meters  and  the  resulting 
grid  would  be  quite  large.  A  continuation  scheme  would  require  far-field 
propagation  Green's  functions  for  a  medium  with  a  bench  on  the  free  surface,  for 
which  there  currently  are  no  analytic  solutions.  Another  advantage  of  the 
reciprocal  approach  is  that  one  calculation  yields  the  excitation  for  many  source 
locations  in  and  near  the  bench  (at  each  grid  point). 

McLaughlin  and  Jih  (1988)  show  that,  if  U(r,s)  is  the  displacement  in  the 
very  far-field  due  to  a  point  explosive  source  along  the  direction  s  at  delayed 
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Figure  1.  Schematic  drawing  of  a  row  of  charges  at  the  face  of  a  quarry  blast. 
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Figure  2.  Schematic  drawing  of  the  reciprocal  theorems  used  in  this  report. 


time  x,  and  if  v(x0,r)  is  the  velocity  due  to  a  plane  wave  incident  at  xc  from 
infinity,  then 

U(t,s)  =  CV  v(x0,r),  (1) 

where  C  is  a  constant  of  proportionality.  Note  that  the  reciprocal  theorem  states 
the  equivalence  between  point  force  sources  and  displacements  in  the  same 
direction  so  that  by  saving  displacements  as  well  as  the  dilatations,  we  can  also 
deduce  the  response  due  to  point  forces  in  the  bench.  A  similar  approach  is  used 
by  Okamoto  (1993).  The  finite  difference  methods,  in  particular  those  for  treating 
the  free  surface  topography,  are  discussed  in  Jih,  et  al  (1988)  and  McLaughlin  and 
Jih  (1988). 

We  have  performed  two-dimensional  finite  difference  calculations  for 
incident  plane  P  and  SV  waves  at  several  angles  of  incidence  as  well  as  for 
Rayleigh  waves  incident  upon  the  bench  from  up-slope  and  down-slope,  as 
shown  in  Figure  3.  The  numerical  grid  is  400  meters  in  the  horizontal  direction 
(with  the  bench  at  the  center)  and  200  meters  in  the  vertical  direction,  with  a 
uniform  grid  spacing  of  2  meters.  The  calculations  were  run  at  500  time  steps  per 
second  for  1001  cycles  (duration  of  0.2  seconds).  The  bench  height  is  18  meters 
and  slopes  at  45°,  approximated  by  9  steps  of  2  meters.  The  calculations  are  valid 
to  about  20  Hz,  above  which  numerical  dispersion  corrupts  the  results.  The  high 
sample  rate  was  dictated  by  the  small  grid  dimension  required  to  define  the 
bench.  For  the  body  wave  calculations,  wave  motion  was  initiated  from  the  right 
side  of  the  grid  along  a  line  inclined  normal  to  the  incidence  angle  passing 
through  a  point  600  meters  from  the  free  surface  in  the  center  of  the  grid.  The 
wave  shape  is  that  if  an  Ohnaka  pulse.  The  Rayleigh  wave  calculations  were 
initiated  500  meters  to  the  left  of  the  bench.  The  initial  motions  had  the  spatial 
and  temporal  properties  of  a  Rayleigh  wave  convolved  with  a  Ricker  wavelet 
(Sheriff  and  Geldart,  1982).  In  addition  to  the  calculations  with  the  bench,  we 
made  comparable  calculations  with  a  flat  surface,  as  references.  Bv  comparing 
the  reference  calculations  to  analytic  solutions,  we  verified  the  usable  bandwidth 
of  the  suite  of  calculations.  Then  by  comparing  the  reference  calculations  with 
those  including  the  bench,  we  isolated  the  effects  of  the  bench. 


2.2.  Results  of  Simulations 
2.2.1  Body  Waves 

At  the  locations  shown  in  Figure  4,  we  monitored  vertical  and  horizontal 
displacement,  dilatation  and  shear  strain  due  to  incident  P  and  SV  waves.  We 
show  in  Figures  5  and  6  the  vertical  displacements  and  dilatations,  respectively, 
due  to  a  P  wave  incident  at  30°  from  the  vertical  traveling  right  to  left.  The 
amplitudes,  indicated  by  both  contour  lines  and  grayscale  intensitv,  are 
normalized  to  the  amplitudes  from  a  reference  calculation  with  a  flat  surface  and 
the  same  incident  P  wave.  The  reference  amplitude  at  a  grid  location  is  the 
amplitude  at  the  same  distance  below  the  free  surface  as  the  location.  Thus,  using 
the  reciprocity  arguments  presented  above,  the  amplitudes  can  be  viewed  as  the 
excitation  at  each  grid  point,  relative  to  the  excitation  in  a  half  space.  The  signals 
have  been  low-pass  filtered  using  an  8-pole  Butterworth  filter  whose  corner 
frequency  is  15  Hz.  As  can  be  seen  from  Figure  5,  the  vertical  displacements 
deviate  very  little  from  1.0  across  the  grid.  The  dilatations  show  much  greater 
variability,  especially  near  the  bench.  Figure  7  zooms  in  on  the  bench,  where  it 
can  be  seen  that  amplitudes  are  reduced  by  as  much  as  a  factor  of  5  on  the  bench, 
and  increased  by  a  factor  of  3  directly  in  front  of  the  bench.  Since  quarry  blasting 
takes  place  behind  the  bench  at  a  horizontal  distance  comparable  to  the  height  of 
the  bench  (Figure  1),  the  effect  of  the  bench  is  to  reduce  the  excitation  of  waves 
from  the  explosive  charges  (the  dilatational  sources)  by  about  a  factor  of  two  for 
waves  traveling  to  the  right  (in  the  sense  in  Figures  6  and  7).  However,  the 
signals  caused  by  movement  of  spall  mass,  represented  by  vertical  and 
horizontal  forces,  are  effected  very  little  by  the  bench,  as  evidenced  by  the  small 
variation  in  displacement  seen  in  Figure  5.  In  Figure  8,  we  show  the  dilatation 
field  for  a  P  wave  incident  at  60°  from  vertical.  The  results  are  similar  to  those  in 
Figure  6,  where  the  incidence  angle  is  30°,  indicating  that  the  excitation  is  not 
sensitive  to  take-off  angle.  The  dilatation  field  for  a  different  corner  frequency,  10 
Hz,  and  for  a  take-off  angle  of  30°  are  shown  in  Figure  9.  Comparison  with 
Figure  7  shows  that  in  the  bandwidth  of  these  calculations,  the  conclusions  are 
not  dependent  on  frequency. 

Note  that  in  Figures  5  through  8  there  is  an  amplification  directly  in  front  of 
the  bench.  Apparently,  the  bench  concentrates  compressive  stresses  at  the  foot  of 
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Figure  7.  Exploded  view  of  the  bench  from  Figure  6:  Dilatations  due  to  a  P 
wave  incident  at  30°  from  the  vertical,  low  pass  filtered  with  comer 
at  15  Hz.  The  dimensions  of  the  plot  are  110  m  x  50  m. 
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Figure  8.  Dilatations  due  to  a  P  wave  incident  at  60°  from  the  vertical,  low  pass  filtered  with  corner  at  15  Hz. 
The  dimensions  of  the  plot  are  140  m  x  50  m. 


Figure  9.  Exploded  view  of  the  bench  showing  dilatations  due  to  a  P  wave 
incident  at  30°  from  the  vertical,  low  pass  filtered  with  comer  at  10 
Hz.  The  dimensions  of  the  plot  are  100  m  x  50  m. 


the  bench  while  reducing  them  behind  it.  The  amplification  is  not  relevant  to 
interpretations  of  observations  of  standard  quarry  blasts  since  it  is  not  usual 
practice  to  set  charges  on  the  quarry  floor.  However,  this  is  a  possible  mechanism 
for  increasing  the  amplitude  of  signals  from  quarry  charges  relative  to  a 
clandestine  bomb. 

We  show  in  Figure  10  the  dilatation  field  for  a  wave  impinging  from  behind 
the  bench  rather  than  from  the  front,  as  in  Figures  5  through  9.  The  incidence 
angle  is  30°.  In  Figure  10,  the  excitation  is  diminished  behind  the  bench,  similar 
to  Figure  7.  The  amplification  at  the  floor  of  the  bench  also  extends  up  the  bench 
slightly. 

The  proximity  of  the  bench  acts  to  produce  SV  waves  from  an  explosive 
source.  Consider  Figures  11  and  12,  in  which  the  dilatation  field  due  to  an 
incident  SV  and  the  rotation  field  due  to  a  P  wave,  respectively,  are  shown.  The 
incidence  angle  is  30°  in  the  figures.  The  strains  in  both  cases  are  the  result  of 
reflections  at  the  free  surface  -  they  are  zero  for  the  incident  wave.  As  with  the 
dilatations  in  Figure  6,  there  is  an  increased  apparent  "source  conversion” 
directly  in  front  of  the  bench  and  a  reduction  behind  the  bench  in  both  Figures  11 
and  12. 

SV  waves  excited  by  a  deviatoric  source  are  enhanced  behind  the  bench,  as 
seen  in  Figure  13,  which  displays  the  shear  strain  field  due  to  an  SV  wave 
incident  at  30°.  Thus,  for  a  source  with  a  deviatoric  component  such  as  a 
compensated  linear  vector  dipole  (CLVD),  the  shear  waves  are  enhanced  behind 
the  bench,  while  the  P  waves  are  reduced. 

2.2.1  Rayleigh  Waves 

Figures  14  and  15  show  the  normalized  dilatation  field  for  Rayleigh  waves 
experiencing  a  step  up  in  topography  and  a  step  down,  respectively.  The  monitor 
stations  were  located  at  the  positions  shown  in  Figure  16.  Again,  by  reciprocity, 
the  values  on  the  plots  can  be  interpreted  as  the  relative  excitation  of  Rayleigh 
waves  due  to  dilatational  sources.  The  Rayleigh  wave  amplitudes  are  normalized 
to  the  incoming  signal  at  the  left  side  of  the  grid.  In  each  case,  going  up-step  and 
down-step,  there  is  a  reduction  of  about  a  factor  of  three  to  five  behind  the  bench. 
The  effect  of  the  step  extends  from  the  step  to  a  distance  about  equal  to  the  step 


14 


Figure  10.  Exploded  view  of  the  bench  showing  dilatations  due  to  a  P  wave 
incident  at  30°  from  the  vertical,  low  pass  filtered  with  corner  at  15 
Hz  (step  down).  The  dimensions  of  the  plot  are  100  m  x  50  m. 
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Figure  12.  Exploded  view  of  the  bench  showing  rotations  due  to  a  P  wave  incident  at  30°  from  the  vertical,  low 
pass  filtered  with  corner  at  15  Hz.  The  dimensions  of  the  plot  are  120  m  x  50  m.  The  amplitudes  can 
be  viewed  as  the  excitation  of  P  waves  by  a  deviatoric  source  at  each  grid  point. 


8 


Figure  13.  Rotations  due  to  a  SV  wave  incident  at  30°  from  the  vertical,  low  pass  filtered  with  corner  at  15  Hz. 

The  dimensions  of  the  plot  are  220  m  x  90  m.  The  amplitudes  can  be  viewed  as  the  excitation  of  SV 
waves  by  a  dilatational  source  at  each  grid  point. 


Figure  14.  Normalized  dilatation  field  for  Rayleigh  waves  experiencing  a  step 
up  in  topography,  low  pass  filtered  with  corner  at  15  Hz.  The 
dimensions  of  the  plot  are  330  m  x  90  m.  The  values  on  the  plots 
can  be  interpreted  as  the  relative  excitation  of  Rayleigh  waves  due 
to  dilatational  sources.  The  Rayleigh  wave  amplitudes  are 
normalized  to  the  incoming  signal  at  the  left  side  of  the  grid. 
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Figure  15. 


Normalized  dilatation  field  for  Rayleigh  waves  experiencing  a  step 
down  in  topography,  low  pass  filtered  with  corner  at  15  Hz.  The 
dimensions  of  the  plot  are  290  m  x  90  m.  The  values  on  the  plots 
can  be  interpreted  as  the  relative  excitation  of  Rayleigh  waves  due 
to  dilatational  sources.  The  Rayleigh  wave  amplitudes  are 
normalized  to  the  incoming  signal  at  the  left  side  of  the  grid. 
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height  for  the  step-down  case  (Figure  15).  This  is  true  for  the  step-up  case  (Figure 
14)  to  the  left  of  the  step,  but  in  addition,  the  amplitudes  are  reduced  to  a  value 
of  0.4  everywhere  on  the  plot  to  the  right  of  the  step.  In  either  case  (step  up  or 
down),  the  source  excitation  near  the  bench  is  significantly  reduced.  The 
corresponding  displacement  fields  are  shown  in  Figures  17  and  18.  The  step- 
down  case  shows  a  reduction  in  amplitude  of  about  20  percent  behind  the  step 
(Figure  18),  while  the  field  is  uniformly  reduced  by  a  factor  of  0.6  to  the  right  of 
the  step  for  the  step-up  case  (Figure  17).  Thus,  Rayleigh  wave  amplitudes  from  a 
vertical  force  near  the  step  will  be  slightly  reduced  near  the  step  if  the  wave 
travels  up-step,  but  will  be  reduced  by  about  40  percent  anywhere  behind  the 
step  if  the  wave  travels  down-step. 
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Figure  18.  Normalized  displacement  field  for  Rayleigh  waves  experiencing  a 
step  down  in  topography,  low  pass  filtered  with  corner  at  15  Hz. 
The  dimensions  of  the  plot  are  320  m  x  90  m.  The  values  on  the 
plots  can  be  interpreted  as  the  relative  excitation  of  Rayleigh  waves 
due  to  vertical  point  forces.  The  Rayleigh  wave  amplitudes  are 
normalized  to  the  incoming  signal  at  the  left  side  of  the  grid. 


3.  Discussion 


Recent  work  by  Stevens,  et  al.  (1993)  has  demonstrated  that  the  long  period 
seismic  waves  radiated  by  an  explosion  near  a  vertical  fre^  surface  mav  be 
significantly  reduced  with  respect  to  an  explosion  under  a  flat  free  surface.  Such 
situations  occur  for  explosions  in  a  mountain  or  ridge.  Stevens,  et  at.  (1993) 
considered  only  the  surface  waves  radiated  by  such  structures.  Here,  we  extend 
the  discussion  to  the  effects  of  a  vertical  free  surface  adjacent  to  an  explosion 
source  such  as  might  occur  in  a  quarry.  Because  of  the  scale  of  typical  quarrv 
faces  (less  than  30  meters)  the  pertinent  wavelengths  of  phases  such  as  P,  Pg,  Lg, 
and  Rg  in  regional  seismic  waveforms  are  much  longer  than  the  quarrv  face.  The 
same  mechanism  that  is  responsible  for  the  reduction  in  surface  wave  excitation 
by  an  explosion  in  a  mountain  applies  to  the  reduced  excitation  of  regional 
phases  by  an  explosion  behind  a  vertical  quarry  face.  Indeed,  we  argue  that  the 
effective  moment  tensor  component  that  describes  the  horizontal  couple 
perpendicular  to  the  face  of  the  quarry  is  reduced.  A  simple  derivation  for  this 
effect  is  given  below.  The  following  mathematical  exposition  provides  a 
theoretical  framework  for  the  interpretation  of  the  simulations  presented  here 
and  a  formalism  for  including  a  source  near  a  bench  in  layered  media 
calculations.  We  show  that  for  wavelengths  longer  than  the  burden  of  the  quarrv 
blast,  we  can  modify  the  affected  component  of  the  explosion  moment  tensor  and 
a  halfspace  Green's  function  may  be  used  for  seismic  propagation  calculations. 
This  approximation  simplifies  the  computation  of  theoretical  seismograms  from 
a  quarry  blast  for  frequencies  of  interest. 

We  are  interested  in  the  effective  radiation  pattern  from  an  explosion  source 
behind  a  quarry  bench.  The  geometry  is  shown  in  Figure  19.  The  point  explosion 
source  is  located  at  rj  i=  0,  a  distance,  h,  to  the  left  of  the  vertical  bench,  h  is 
referred  to  as  the  burden  of  the  quarry  blast.  Typical  quarry  blast  practice  is  to 
choose  a  burden  between  30  percent  and  100  percent  the  height  of  the  bench. 

The  seismic  displacement  at  a  location,  X,  from  an  explosion  source  with 
moment  tensor,  Mjj  =  8^  M0  at  location  f)  may  be  written  as 

U,(x)  =  ^ckkpq(;7)G1[).q(X,>j),  (1) 
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b  =  Bench  Height 

h  =  Burden  or  Source  Distance  to  Bench  Face 


Figure  19.  Bench  geometry. 
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« 

where  Uj^Xj  is  the  i'th  component  of  displacement  at  location,  X,K  =  A  +  2q  /  3 
is  the  bulk  modulus  at  the  source  location,  C^pq(  77)  is  the  elastic  tensor,  and 

Gip,q(X  77)  is  the  q’th  derivative  of  the  Green's  tensor  response  for  the 

i'th 

• 

component  of  motion  at  X  given  a  force  in  the  p'th  direction  at  location  r) . 

Note  that  by  reciprocity. 

• 

Ckkpq(  ^)Gip,q(*,  *l)  =  CTli(  V)  +  G-nk  n)  +  ^33(  *?), 

(2) 

where  Oi  1, 022  and  033  are  stress  terms  at  location  77  due  to  a  displacement  in  the 

• 

i'th  direction  at  location  X.  We  now  make  use  of  this  reciprocal  problem,  where 

Oij  is  the  stress  tensor  at  i)  due  to  a  displacement  at  X.  Similarly,  we  will  write  u^ 

as  the  displacement  at  f\  due  to  a  displacement  at  X.  Dependence  on  the  location 

X  and  the  i'th  direction  of  motion  is  suppressed  for  clarity  of  notation  in  the 

• 

following  discussion. 

Now,  we  know  on  the  vertical  bench  (a  vertical  free  surface)  that  on  =0,  and 

we  expand  the  131  dependence  of  On  in  a  Taylor’s  expansion  from  the  bench. 

• 

~  -  o'il(7?i  -  h)  h  |i?i=h  -  h  ^  |r?1=h 

(3) 

We  have  the  equation  of  motion  which  states  that 

# 

*,>=h+  *,'>=h=pUl(h) 

(4) 

together  with  the  free  surface  boundary  conditions  on  the  bench  (at  r/i  =  h  ) 

# 

(5) 

hence  we  can  state  for  the  horizontal  equation  of  motion  on  the  bench 

» 

7nUhi=h  =PUi(7i  =h) 

rJTh 

(6) 

• 
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Therefore, 


^n(»?i=0)  =  -hpui(n1=h)  (7) 

We  approximate  iii(oi  =  h)  by  the  halfspace  solution,  u/HS^(r?]  =  h),  and  we  can 
write 

On(*?i  =0)  =  -hpu/HS)(ni  =  h) .  (8) 


This  says  that  the  an  stress  from  a  displacement  source  at  location X  is 
proportional  to  the  distance  from  the  bench,  h,  and  density,  p. 


How  much  is  the  explosion  behind  a  bench  reduced  with  respect  to  an 
explosion  in  a  halfspace?  We  can  estimate  this  effect  by  writing  the  expression  for 
the  halfspace  an  component  of  stress  in  terms  of  the  halfspace  displacements. 


oil 


,  WHS)  / 


<WHS)  |  <9u3(hs) 
dr]2  dTj3 


(9) 


To  estimate  the  ratio  F  =  -^jgjwe  neglect  the  second  term  to  get 


'll 


on 


(HS) 


pa‘ 


<WHS) 

dm 


(10) 


We  note  that  for  plane  P  waves  propagating  in  the  rji  direction  we  can  make  the 
approximation  that 


dm 


(ii) 


so  we  have  that 


r  = 


an  _  phta2u/H^  _  wh 


On 


(HS) 


paaj Ui 


(HS) 


=  kh 


a 


(12) 


For,  h  =  10  m,  a  =  2000  m/s,  we  have  T(lOHz)  =  0.31,  T(5Hz)  =  0.15,  and  T(lHz)  = 
0.031.  This  would  indicate  that  the  P-waves  radiated  perpendicular  to  the  bench 
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will  be  substantially  reduced  with  respect  to  the  halfspace  for  much  of  the 
bandwidth  of  regional  seismograms. 

However,  it  is  clear  that  the  radiation  pattern  from  an  explosion  behind  the 
quarry  bench  will  depend  upon  the  slowness  and  the  azimuth  of  the  radiated 
wave.  That  is  to  say  that  T  =  T^k,hj.  A  complete  treatment  of  this  problem 

therefore  requires  the  solution  of  the  complete  scattering  problem  using 
numerical  techniques  as  in  the  previous  section. 

To  summarize,  for  kh  «  1  the  Mu  component  of  an  explosion  behind  the 
bench  can  be  modified  and  halfspace  Green's  functions  can  be  used  to  propagate 
the  seismic  waves.  To  formalize  this,  we  re-write  Equation  (1)  as, 

Ui(x)  =  H(rn(k,h)ckknGn,i(HS)(x,'j)+  .  <i3) 

r22(k,h)Ckk22Gi2,2  <HS>(x,  ri)+ r33(k,h)ckk33Gi3,3(HS)(x,  77)) 

For  kh  «  1,  we  have  argued  that  Tn  (k,h)  *  kh  «  1,  while  T^k,!^  =  1  and 

r33^k,hj  ~  1.  And  we  approximate  Equation  (13)  by, 

ui(x)=||-(khckklIGiu(HS)(X,f;)  (14) 

+  Ckk22  Gi2,2^HS^(X'  Vj  +  Ckk33  Gi3,3^(X, 
which  may  be  more  simply  written  for  an  isotropic  medium, 

U4(X)  -  Mc(  khGiu(HS>(x,  77}  +  Gi2,2(HS)(x, 77)  +G13,3(HS)(X,  tj))  •  (15) 

It  should  be  readily  apparent  from  this  expression  that  the  radiated  seismic 
waves  will  have  an  azimuthal  radiation  pattern  and  the  source  will  radiate 
S  H  waves. 
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4.  Conclusions 


We  have  computed  the  strain  and  displacement  fields  due  to  plane  waves 
incident  on  a  uniform  solid  with  a  topographic  bench  at  the  surface.  Using  the 
reciprocity  arguments  discussed  above,  this  allows  us  to  map  the  relative 
excitation  of  waves  due  to  point  dilatational  and  force  sources  at  and  near  the 
bench.  We  find  that  the  excitation  of  waves  on  the  bench  where  quarry  blasting 
would  occur  are  influenced  relative  to  the  excitation  of  waves  in  a  flat  half-space, 
in  the  following  ways. 

1.  Seismic  signals  from  point  explosion  (dilatational)  sources  located  behind  the 
quarry  face  can  be  reduced  substantially  in  amplitude,  while  seismic  waves 
from  sources  on  the  quarry  floor  may  be  amplified.  The  source  amplitude 
variations  occur  for  sources  within  about  one  bench  height  of  the  quarry  face. 
The  numerical  simulations  presented  in  this  study  can  be  used  to  estimate 
and  correct  for  this  effect. 

2.  An  extended  source  in  the  standard  configuration  (charges  placed  a  distance 
of  1/2  the  bench  height  from  the  quarry  face  and  extending  to  the  depth  of 
the  quarry  floor)  averages  these  variations  leading  to  a  source  that  is  reduced 
slightly  in  amplitude  relative  to  a  source  in  a  half  space. 

3.  The  amplitude  variations  occur  because  a  dilatational  source  such  as  an 
explosion  is  strongly  affected  by  the  presence  of  a  free  surface  close  to  the 
source.  Seismic  signals  from  point  forces  are  affected  very  little  by  proximity 
to  the  quarry  face. 

4.  These  results  are  nearly  independent  of  the  takeoff  angle  of  the  seismic  signal. 

5.  Rayleigh  wave  amplitudes  are  reduced  by  up  to  40  percent. 

In  terms  of  modeling  quarry  blast  ground  motions,  such  as  in  Barker,  et  al.  (1993), 
this  study  suggests  that  when  computing  regional  phases  such  as  Pn,  Pg,  Lg 
using  plane  layered  models,  the  Mn  component  from  the  explosive  charges 
should  be  reduced  by  a  factor  of  2  or  3  while  the  forces  from  the  spall  phases 
need  not  be  scaled.  The  contributions  to  regional  phases  from  spall  were  found  in 
Barker,  et  al.  (1993)  to  be  comparable  to  those  from  the  explosive  charges.  The 
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results  here  argue  that  the  explosive  component  of  the  blast  be  reduced  up  to  a 
level  about  50  percent  that  of  the  spall  component.  Therefore,  in  the  far-field,  for 
wavelengths  longer  then  the  quarry  bench  height,  the  spall  components  of  a 
quarry  blast  are  expected  to  exceed  the  explosive  part  of  the  quarry  blast  for 
wave  radiated  perpendicular  to  the  quarry  bench. 

Furthermore,  the  apparent  reduction  of  Mu  relative  to  M22  and  M33  implies 
an  apparent  deviatoric  source  that  will  radiate  SH  waves  from  an  explosion 
behind  the  bench.  The  enhanced  SH  and  SV  radiation  contribute  to  the  Lg 
excitation  by  quarry  blasts  relative  to  normally  buried  and  overburied  explosions 
that  are  not  adjacent  to  vertical  free-surfaces  such  as  a  quarry  face. 
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